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Abstract 

A condition for determining the holonomy group associated to a 
principal bundle with connection is proven. Within the context of 
Holonomic Quantum Computation, this condition aids in determining 
the universality of the system. We apply the theory to the product 
bundle arising from the full two-qubit system of holonomic compu- 
tation with squeezed coherent states and conclude that the control 
algebra generates the Lie algebra u(4) . 

1 Introduction 

Geometric phases have long been a source of fascination and insight into 
classical and quantum physical theories |l| . In recent years, they have proven 
to be useful in describing the dynamics and control of certain nonholonomic 
mechanical systems with symmetry. Inspired by the appearance of geometric 
phases in biology, engineers have sought to create motion in machines via 
cyclic variations in shape space. These endeavors and the characterization 
of optimal trajectories [[3], |J remain active areas of research. 

Most recently, the geometric phase has been realized as a way of construct- 
ing logic gates in a quantum computer ||. Holonomic Quantum Computa- 
tion (HQC) employs non-Abelian geometric phases (holonomies) for the pur- 
pose of quantum information processing. This paper is an elaboration of . 
There it was reported that the Optical Holonomic Computer 0, ||, |9], |10[ [11 



is universal by showing that the holonomy group of the two-qubit interaction 
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model contains SU(2) x 577(2). These transformations act non-trivially on 
the standard computational basis and support, in particular, SWAP gates. 
Coupled with the ability to perform single qubit rotations, the claim of uni- 
versality follows by appealing to standard results from quantum computation 
P. 

These notes and completely solve the problem posed in the recent 
work of K. Fujii |9|, [H| |IT| . That is, for the two-qubit system, we explic- 
itly construct a set of holonomic transformations that span the Lie algebra 
w(4) and thus attain universality of the system. Note however that we do 



not carry out the program proposed in |TTJ. In Pachos-Chountasis [[7] and 
Pachos-Zanardi || , the universality of the model was shown by constructing 
a universal two-qubit logic gate. This £7(4) transformation was integrated di- 
rectly by choosing a loop in a 2-submanifold of parameter space along which 
the connection components commute but with non-vanishing curvature. For 
these restricted loops, the non-Abelian nature of the path ordered integral 
may be avoided. However, that transformation was not an element of the 
holonomy group as calculated by the application of the Ambrose- Singer The- 
orem in ||10|| . 

In ||, we stated simple conditions for determining the holonomy group of 
a principal bundle with connection. These conditions are seen to be a corol- 
lary to Chow's theorem from control theory. The obvious relevance to HQC 
is that the holonomy group represents the set of controlled transformations 
attainable for manipulating quantum information. 

In this section, we review some geometric constructions and provide some 
details regarding the corollary stated in ||. This proof is not new. A similar 
proof can be found in Q and, while not explicitly stated, all the necessary 
ingredients can be found in [13|j. 



1.1 Principal and associated bundles 

Let Q be a principal G— bundle defined by the free (right) Lie group action 
$ : Q x G — > Q and the canonical projection Q — > Q/G = M. If U is a 
neighborhood of M, then Q is locally diffeomorphic to the product U x G. A 
smooth map a : U — >Q, such that 7r o a = idu is called a local section over 
U. 

Given any principal bundle Q /G and an action of the Lie group G 
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on a vector space V, we may define 



p: QxV 
(?,«) 



V 



(1) 



where q£Q,v£V,g£G and p is a representation of the group G. 
The associated vector bundle E p with typical fiber V is the set of equivalence 
classes under this action 



where points (q,v) and (qg, p~ 1 (g)v) are identified. Let [(q, v)] denote a 
representative of the equivalence class. 

In the mathematical theory of HQC, there are two associated bundles 
of interest. First, the quantum code is realized as the C n vector bundle 
associated to the pullback of the universal bundle by the classifying map 
H [HJ . Second, the adjoint bundle defined by taking V = g (the Lie algebra 
of G) and p as the adjoint action of a Lie group on its Lie algebra is important 
in deriving conditions for universal computation as we prove in Section 1.3. 

We now state an important correspondence concerning the associated 
adjoint bundle. First, recall that a Lie algebra valued differential form uj is 
74d-equivariant if 



If uj is a Lie algebra valued function f on Q, then Ad-equivariance is written 
as f(qg) — Ad(g _1 ) f '(g) , g G G. The sections s of the adjoint bundle 
associated to a principal bundle Q are in one-to-one correspondence with 
Ad-equivariant functions on Q. To any Ad-equivariant function / we may 
define a local section s/ by s/(p) = [(g, /(g))], where q is any point in the 
fiber over p. This definition is independent of the choice of g because of Ad- 
equivariance. Conversely, given a local section s : U ->Q x Ad g, we define a 
map f s : Tr -1 ^) such that s(p) = [(pj s (p))] 0- 

1.2 Connections and covariant derivatives 

Recall that for any £ 6 q, the R-action $ e xp(t€) Q defines a curve though q e Q. 
The infinitesimal generator of the action corresponding to £ is defined as the 



E p = Qx p V = (Qx V)/G, 



(2) 




(3) 
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vector field 



= T+ $ex p(t£) 1- 



(4) 



t=o 



The vertical subspace V q Q C T q Q is defined to be the subspace of T q Q that 
is tangent to the fiber over p — 7r(g) , G p . By the previous definition, we have 
the identification V q Q = q. 

A connection on a principal G-bundle A on Q is an Ad-equivariant Lie 
algebra-valued one-form A : TQ — > q such that A(£ q ) = £. The horizon- 
tal space H q Q is the linear space H q Q = {X G T q Q\A q (X) = 0}. These 
definitions provide the splitting 



of the tangent vectors into horizontal and vertical components. Note that 



The map T q n : H q — > T^t q \M is an isomorphism from the horizontal space 
to the tangent space of the base space. Thus a curve q{t) e Q defines a 
curve in the base by specifying a tangent vector at each point it(q(t)) = pit). 
The properties of the connection and the uniqueness theory of ODE's provide 
the reverse procedure of reconstructing a unique curve in the horizontal space 
given a curve in the base called the horizontal lift |15| . Denote the horizontal 
lift of the vector X e T p M where p = ir(q) by X^ and hX q as the horizontal 
component of a vector in T q Q. 

Let a be a local section over a neighborhood U of M that contains a 
curve 7. Then the horizontal lift 7 can be expressed as 7(t) = cr (7 (t)) <?(£), 
where = g( , ~f(t)) G G. The local connection one-form A is defined as the 
pullback A = a* (A). The local connection may be written in terms of group 
elements by the equation 

*f = -AM. (?) 

Choosing the initial condition, g{0) = idc, the solution may be obtained by 
the path ordered integral 



T q Q = V q Q © H q Q 



(5) 



V q Q = KerTgTT, and H q Q = Ker^4 9 . 



(6) 
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(8) 
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Given a connection on a principal bundle, we may define the covariant 
derivative of sections of associated bundles. Let E g = Q x Ad g be the adjoint 
bundle associated to Q as defined above. Let 7 be a base curve with tangent 
vector X at 7(0). If {T a } is a basis of g, then a section s of E s may be 
obtained by forming s = [(a(p),T(p))], where T{p) = f a {p)T a and a is a 
local section over a neighborhood U of p. The covariant derivative Dxs of s 
along X G TM is defined by [[15 



Dxs = [(7(0) , jAd g - Ht) T(t)] J] 



t=o 



The previous definitions imply, 

= (Wo) 



[(a(0)g(0) , |Ad ff - 1(t) T(t)] t=Q , 
Ad^Ad^T^J] 
[K0),f + [A(X),T(*)]| )]. 



(9) 

(10) 

(11) 
(12) 



Covariant derivatives of sections of the adjoint bundle provide a way of 
computing Lie derivatives of Ad-equivariant functions on the total space by 
the correspondence 



D x s f = X h f 



(13) 



where Sf is the section of the adjoint bundle corresponding to the Ad- 
equivariant function / [ p~3]j . 



1.3 Infinitesimal holonomy 

The horizontal lift of any closed curve in the base space with initial point p 
maps the fiber over p to itself and thus corresponds to the group translation 
qg, q G 7r _1 (p) and g G G. In this way, we can associate to closed curves 
in the base space elements of the structure group G. The set of all such 
group elements forms a subgroup of G called the holonomy group. With the 
restriction to sufficiently small loops homotopic to zero, we may define the 
infinitesimal holonomy group. The infinitesimal holonomy group is in general 
a subgroup of the full holonomy group, however for analytic connections over 



connected base manifolds, these notions coincide 13 
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The holonomy group is determined from the curvature. The curvature T 
associated to A is by definition the Lie algebra-valued two-form on Q such 
that 

F{X 1 ,X 2 ) = dA(hX 1 ,h X 2 ) (14) 

where d denotes exterior differentiation. Using the identity 

dA(X 1: X 2 ) = X 1 A(X 2 ) - X 2 A{X l ) - A([Xx ,X 2 \) (15) 

and the fact that the connection form evaluates to zero on horizontal vectors, 
we get the relation !F(Xi, X 2 ) = — A([h Xi , h X 2 ]). Roughly, the curvature 
measures the non-commutativity in the total space of horizontally lifted vec- 
tor fields defined in the base space . Note that T is Ad-equivariant since A 
is Ad-equivariant and pullback and exterior differentiation commute. 

The Lie algebra corresponding to the infinitesimal holonomy group is 
spanned by the g-valued Ad-equivariant function / given by 

/ = x l x t- x --- x in^i^i) as) 



where Xj G TM and 2 < k < oo |16|. By the correspondence (|T3D, the 
pullback through a local section a of this function may written in terms of 
covariant derivatives as 

D Xik D Xiki ...D Xi3 F(X i2 ,X h ) (17) 

where F is the local form of the curvature defined with respect to a local 
section F = a* (J 7 ). 



2 HQC with squeezed coherent states 

In this section, we revisit the mathematical framework of Optical Holonomic 
Computation. See [0, ||, ||, [1(J for the many details omitted here. We refer 
to this setup as HQC with squeezed coherent states since the properties of 
interest in quantum computation, e.g. universality, entanglement, etc., rely 
only on the squeezing and displacing operators and not directly on the degen- 
erate Hamiltonian used to encode the qubits. Therefore, the mathematics 
describing this computing scheme will be shared by any quantum system 
whose states can be displaced and squeezed. See for example the remarks 
concerning trapped ions in ||. 
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Let Hi = {\ati) , a G N} be a Fock space with number operator n, . 
The Hamiltonian H l = nj(nj — 1) encodes the z-th qubit in the degenerate 
subspace {|0j) , |lj}} . Single mode squeezing and displacing operators control 
the single qubit 

S l (/ii) = exp (/lidf - pid?) D l {\i) = exp (Aid] - A^) (18) 

where \ii , \ G C and a« and a\ are the annihilation and creation opera- 
tors respectively. To obtain controllable interactions among product states 
we employ the Hamiltonian H 12 = n\{n\ — 1) + n 2 (n 2 — 1) with degener- 
ate subspace {|00), 1 01) , 1 10) , |11)} and two mode squeezing and displacing 
operators 

M(C) = exp (Cal4 - Coia 2 ) N(£) = exp (£a{a 2 - fonaj) (19) 

where f , £ G C. Set U 2 = D(X)S(ji) and U A = N(£)M((). 

The geometry of the U(n) holonomy problem [§, [T^O]] is briefly described 
as follows. Define the manifolds 

st n (H) = {v = ( Vl , . . . , v n ) g n x . . . x n | v*v = u nxn } (20) 

Gr n (H) = {X G B(H) I X 2 = X , X f = X , trX = n} . (21) 

where B(7i) denotes the set of bounded linear operators on TC. These mani- 
folds are known as the Stiefel and Grassmann manifolds respectively. They 
form a principal bundle with the (right) U{n) action on St n (W) and the 
projection n : St n (H) — > Gr n (H) given by 7r(V) = VVK Denote this U{n)- 
bundle by P n . Let M C C 2 be a parameter space and define the classifying 
map II n : M — > Gr n (H) given by 

n 

n n ( v )=u n ( v )(J2 v j v ])ul(v) (22) 

j=i 

Now form the pullback bundle Q n = II*P n . The quantum code is then re- 
alized as the C n - vector bundle associated to Q n || [L0|, [15[]. To ease the 



notation, we consider only the manifolds P n and Q n . So then, the full two- 
qubit system consisting of single qubit rotations on each qubit and interac- 
tions among them is modelled by the pullback bundles corresponding to the 
product bundle 

St 2 (W0 x St 4 (Hi®ft 2 ) x St 2 (W 2 ) 

1 [ W 

Gr 2 (Wi) x (Gr 4 ) int (^i®H 2 ) x Gr 2 (ft 2 ) 
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The local connection one-form is given by the Wilczek-Zee connection [17 
with matrix elements 

A^^(g\U^ V )-^-U( V )\g). (24) 

The connection components are listed in Appendix A with respect to the 
labelling of parameters Ai = Xi + iy\ , \i\ = rie tdl , A 2 = x 2 + iy<i , = 
r±e l6i , ( = r 2 e ie2 and£ = r 3 e l8:i . 

2.1 Reduction 

The local curvature components are listed in Appendix B. The holonomy 
group of the two-qubit interaction model is SU{2) x SU(2) xU(l) ||. These 
transformations leave invariant the subspaces {[00) , |11)} and {[01) , |10)} 
in Tii Cg) Ti,2 . By the reduction theorem for connections ]n|, the connection 
is reducible to a su{2) x su{2) x w(l)-valued connection and we may reduce 
the total space to St2,4(7ii <8> TC 2 ) ■ To determine the reduced base manifold, 
we form the quotient 

U{i) ~ SUi4) -SGr^CW. (25) 



SU(2) x SU(2) x U(l) SU{2) x SU(2) 

In a similar manner, we can reduce the bundles St2(7ij) 1— > Gr2(7ij) cor- 
responding to the single qubit rotations. The U(2) holonomies act in the 
product space Hi ® H2 as U(2) <8> 1 and 1 <8> U{2). The bundles reduce to 
St2,4(Wi ®W 2 ) ^ Gr 2i 4(7ii ®H. 2 )- For the full two-qubit system, we have the 
reduced product bundle 

St 2 ,4Cftl®ft 2 ) X St 2 , 4 (ftl <g>ft 2 ) X St 2 ,4(Hl®^2) 

1 l (26) 

Gr 2j4 (Wi ® W 2 ) x (Gr 2 , 4 ) int (^i®^ 2 ) x Gr 2 , 4 (Wi ® H 2 ) 

This structure is similar to the product bundle for the generic two-qubit 
system proposed by Pachos fl8fl. 



2.2 Control algebra 

The standard argument to prove universality of a quantum system is to 
demonstrate the ability to perform all single qubit rotations on each qubit 
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and a non-trivial two-qubit transformation. The implicit argument is then 
that the control algebra formed by these matrices and their Lie brackets 
generates the Lie algebra su(n) or u(n) . In this section, we carry out this 
argument for HQC with squeezed coherent states. The matrices we refer to 
can be found in Appendices B and C. 

The holonomy algebra for each single qubit system is u{2) . 

u{2) <g> 1 : span R {F xm <g> 1 , F ri6l <g> 1 , F yiri 1 , F Xiri <8> 1} (27) 
1 <g> u{2) : span R {l <g> F X2y2 , 1 <g> F r4 e 4 , 1 <g> F y2J . 4 , 1 <g> F X2r4 } (28) 

The holonomy algebra for two qubit interactions system is 
su(2) x su(2) x u(l) . 

span R {F r2r3 ,F r2 e 2 , F r2 g 3 , F r303 ,Da_F r2d2 ,DaF r , 2 g 2 ,DdDaF r202 } 

dr2 dQ'2 d&2 

(29) 

From these matrices, we choose the subset linearly independent over the reals 

Ci = {F ridl <g> 1 , F yiri <g> 1 , F xm <g> 1} (30) 
C 2 = {1 <g> F X22/2 , 1 <g> F OTr4 , 1 <g> F X2r J (31) 
C12 = {F rar3 ,F r . 2 o 2 ,F r . 2 g 3 ,D_a_F r2 e 2 , D_a_D_a_F r2d2 } (32) 

^^2 ^2 

Additional independent matrices can be obtained by taking Lie brackets 
from the set {C\ U C 2 U Ci 2 } . In principle, these transformations can 
be constructed in the usual infinitesimal way. For example, to construct 
[F xm <8> 1 , F r2T3 ] traverse a small loop in the plane (xi , ri), then traverse a 
loop in the plane (r 2 ,r 3 ), then traverse the loop backwards in (xi , ri), and 
finally traverse the loop backwards in (r 2 , r 3 ) . One choice is given by 

{ [ F xiri <S> 1 , F r2r3 J , [ 1 <S> F X2Ti , F r2r3 J , [ F yiri ® 1 , F r2 Q 3 J , 

[ 1 <g) F y2r4 , F r26 / 3 ] , 1 ® F^^ , [F xin <S> 1 , F r2T3 \ } (33) 

Note that these matrices are independent with each other and the set 
{Ci U C 2 U Ci 2 } over the reals. 

The 16 matrices C = {Ci U C 2 U C i2 U C^} span the Lie algebra «(4) . 
Thus, at least in principle, we have obtained complete controllability over 
the two-qubit system. 
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Appendix 



A Connection components 
A.l Single qubit 



A Xl — 



— iyi —(cosh 2ri — e tdl sinh 2r\) 

cosh 2ri — e^ 6 * 1 sinh 2ri — «/i 

ixi « (cosh 2ri + e ldl sinh 2ri) 

i(cosh 2ri + e"^ 1 sinh 2ri) ixi 



0i 



1 
3 



-(cosh 4?"! — 1) A ri 








A 



—iyi —(cosh 2r4 — e ldi sinh 2r 4 ) 

cosh 2r4 — e - * 6 * 4 sinh 2r4 — iyi 

1X2 i(cosh 2r 2 + e l9i sinh 2r 4 ) 

i(cosh 2r 4 + e - * 6 * 4 sinh 2r 4 ) ix 2 

3 7( cosh4r 4-l) A 4 = go 



A.2 Two-qubit 



A r2 — 









e i(>2 



— e~ 









A r3 — 







e Wz 





— e 








(2 cosh 2 r 2 



^0 2 = 



103 



" 





e 


-102" 


























f 


2 








"0 








0" 








e -i0 3 








e m 





















- sinh 2r 2 + 



cosh 2r 2 sin 2r 



10 

2 

2 

3 



1 







-(cosh 2r 2 



o 
o 



i sm' 
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Holonomy algebra 



.1 u{2) 



xiyi 



F = 

1 xiyi 



y\r\ 



+ 



+ 




1 

-1 



-1 

-1 


0' 

1 



F ri9i - 



1 
3 



% sinh 4ri 



2i(e iei cosh 2r l + sinh 2r ± ) 

2i (e~ idl cosh 2r\ + sinh 2r x ) 
2(e idl cosh 2ri - sinh 2r\) 
2(e- i01 cosh2rx - sinl^rx) 



1 = 







10 

1 



Ai 



yiri 













>101 



1 



1 = 



10 

10 

3 

3 



i sinh 4ri 








2i(e iei cosh 2r x + sinh 2r ± ) 



+ 



1 = 









-10 

-1 



-1 










2i(e ' Wl cosh2ri + sinh2ri) 








+ 







10 

10 



2(e idl cosh2n - sinh2ri) 



2{e~ idl cosh2rx - sinh2r!) 
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B.2 1(g) u{2) 



X2V2 



Ai 



p — 

1 V2T4 



+ 



P — 



+ 




1 

-1 



-1 

-1 


0' 

1 



r 4 4 



1 
3 



% sinh 4r4 



2i(e if?4 cosh 2r 4 + sinh2r 4 ) 

2i (e~ m cosh 2r 4 + sinh 2r 4 ) 
2(e idA cosh2r 4 - sinh2r 4 ) 
2 (e~ m cosh 2r 4 - sinh 2r 4 ) 



X2D2 



10 





1 



T404 



10 

3 

10 

3 



i sinh 4r 4 



V2T4 



+ 



p — 

1 X2T4 



+ 



0-100 



-1 





-10 


0-10 



2i(e i64 cosh 2r 4 + sinh 2r 4 ) 



2i (e i04 cosh 2r 4 + sinh 2r 4 ) 



0-100 



-1 



0' 

10 



10 



2(e if?4 cosh 2r 4 - sinh 2r 4 ) 



2 (e~ idA cosh 2r 4 - sinh 2r 4 ) 
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B.3 su(2) x su{2) x u{l) 



F = 



"0 














_ e -i03 







































e i0 3 















2 sinh 2r 2 F r202 = 





10 

10 

2 



2i sinh2ro 



i sin 2r 3 sinh 2r 2 F rs g 3 


















-1 














1 


















% sin 2r3 sinh 2 2r 2 , 



D a R 



Dd_F r202 

dr 2 









J8 2 



T2&2 ~ 















je 2 



















-102 









4i sinh 2r 2 



2 sinh 2 2r 2 





10 

10 

2 



4i cosh 2r 2 , 



D a D a F rn f> n 

ae 2 ae 2 



1 











-1 



2isinh 3 r 2 + 



' e" if?2 ' 





e td2 



2i sinh 2r 2 cosh 2r 2 
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C Higher order brackets 



+ 



[ 1 ® F X2r4 , F r2V3 J 



+ 



[F yiri ® 1 , F r2 e 3 ] 



+ 



[ 1 <g> Fj, 2r4 , F r2 3 ] 



+ 



0-100 



1 



0' 

10 


0-10 

10 

-1 









-10 

10 

10 



-1 



0' 

-10 



10 

10 

-1 





0' 



-10 

10 



4^(01+03) sinh 2r 2 ( cosh 2n - e~ iei sinh 2n) 



4e -i(ei+0 3 ) sinh 2r2 ( cos h 2 ri - e iei sinh 2rJ 



4e -i(e 3 -0 4 ) sinh 2r2 ( cos h 2r 4 - e ~ if?4 sinh 2r 4 ) 



4e i(0 3 -04) sinh 2r 2 ( cosh 2r 4 - e if?4 sinh 2r 4 ) 



2e m+e s ) sin 2rs sinh 2r 2 ( cosh 2n + e~ iBl sinh 2n) 



2 e -*(0i+03) sin 2r 3 sinh 2r 2 ( cosh 2r 1 + e if?1 sinh 2r : ) 



2e -m-d4) gin 2r3 sin h 2r2 ( cosh 2r 4 + e~ id4 sinh 2r 4 ) 



2e i(0 3 -04) sin 2rg sinh 2r2 ( cosh 2r 4 + e i04 sinh 2r 4 ) 



14 



-10 
10 
10 
-1 



5 = 16i sinh 2r 2 ^ — sinh 2r x ( cos(0 3 — 6* 4 ) cosh 2r 4 + cos 9 3 sinh 2r 4 ) 
+ cosh 2ri ( cos(6 l 1 + #3 — 6> 4 ) cosh 2r 4 + cos(#i + 6> 3 ) sinh 2r 4 ) j 
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